Preprint typeset using I4TgX style emulateapj v. 10/09/06 



HDE 245059: A WEAK-LINED T TAURI BINARY REVEALED BY CHANDRA & KECK 



C. Baldovin Saavedra 12 , M. Audard 1,2 , G. Duchene j ' 4 , M. Gudel 3 , S.L. Skinner , F. B. S. Paerels', A. Ghez 8 , and C. 

McCabe 5 



,3,4 



ON 
O 

o 

(N 

X> 
<D 
IX, 

o 

(N 
P? 

43 
6 



> 

m 
in 
m 

(N 
O 
On 
O 



13 



ABSTRACT 

We present the Chandra High Energy Transmission Grating Spectrometer (HETGS) and Keck observations 
of HDE 245059, a young weak-lined T Tauri star (WTTS), member of the pre-main sequence group in the 
A Orionis Cluster. Our high spatial resolution, near-infrared observations with Keck reveal that HDE 245059 
is in fact a binary separated by 0."87, probably composed of two WTTS based on their color indices. Based 
on this new information we have obtained an estimate of the masses of the binary components; w 3M Q and 
w 2.5M for the north and south components, respectively. We have also estimated the age of the system to 
be » 2 - 3 Myr. We detect both components of the binary in the zeroth order Chandra image and in the grating 
spectra. The lightcurves show X-ray variability of both sources and in particular a flaring event in the weaker 
southern component. The spectra of both stars show similar features: a combination of cool and hot plasma as 
demonstrated by several iron lines from Fe XVII to Fe XXV and a strong bremsstrahlung continuum at short 
wavelengths. We have fitted the combined grating and zeroth order spectrum (considering the contribution of 
both stars) in XSPEC. The coronal abundances and emission measure distribution for the binary have been 
obtained using different methods, including a continuous emission measure distribution and a multi-T approxi- 
mation. In all cases we have found that the emission is dominated by a plasma between ~ 8 and ~ 15 MK a soft 
component at ~ 4 MK and a hard component at ~ 50 MK are also detected. The value of the hydrogen column 
density was low, Nh ~8x 10 19 cm" 2 , likely due to the clearing of the inner region of the A Orionis cloud, 
where HDE 245059 is located. The abundance pattern shows an inverse First Ionization Potential (FIP) effect 
for all elements from O to Fe, the only exception being Ca. To obtain the properties of the binary components, 
a 3-T model was fitted to the individual zeroth order spectra using the abundances derived for the binary. We 
have also obtained several lines fluxes from the grating spectra. The fits to the triplets show no evidence of high 
densities. In conclusion, the X-ray properties of the weak-lined T Tau binary HDE 245059 are similar to those 
generally observed in other weak-lined T Tau stars. Although its accretion history may have been affected by 
the clearing of the interstellar material around A Ori, its coronal properties appears not to have been strongly 
modified. 

Subject headings: stars: abundances — stars: coronae — stars: pre-main sequence — stars: individual (HDE 
245059) — X-rays: stars 



1. INTRODUCTION 

T Tauri stars (TTS) are optically revealed pre-main se- 
quence stars (PMS) of low mass (M» ~ 0.2-3 M Q ), whose 
interiors are fully convective and powered principally by grav- 
itational contraction rather than by nuclear reactions. Low- 
mass PMS stars are classified optically in two types: classical 
T Tauri stars (CTTS) and weak-lined T Tauri stars (WTTS). 
The CTTS present strong Ha emission lines, a signature of 
accretion, and infrared excess, revealing a dusty circumstellar 
disk. On the other hand, WTTS present weaker Ha emis- 
sion lines and little or no infrared excess, which is generally 
interpreted as a sign that active accretion has significantly de- 
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creased and the disk has become optically thin. PMS stars are 
magnetically active and rotate slowly during the early stages 
of their e volution and spin up as the y contract to the main 
sequence ( Stau ffer & Hartmann 1119861) . Recent evidence sug- 
gests that there is also a population of low-mass PMS stars 
(0.3 < M/Mq < 1) that seem to contract towards the zero-age 
main sequence (ZAMS) at constant angular velocity during 
the first 3-5 M yr after they begin th eir evolution down the con- 
vective tracks (Rebu lTet al. 1 2004). indicating there should be 
a mechanism extracting angular momentum in these PMS ob- 
jects. In CTTS the presence of disk and magnetic fields is be- 
lieved to keep the rotation rate low, while the short rotational 
period in WTTS may be due to the dispersal of the circum- 
stellar disk. But observational evidence has been sometimes 
confusing and the relation betw een accretion an d angular mo- 
mentum loss cannot be assured (lBouvierll2007l) . 

Both type of stars are strong X-rays emitters that 
were first detected w it h the Einstein satellite, e.g. 
Feigel son & De Campli I (119811) and subsequent obser- 
vations with ROSAT (e.g. Feigelson et al. 1993, Neuhauser 
et al. 1995). Spectral measurements of X-rays from YSOs 
reveal emission from continuum and lines from an optically 
thin plasma. Most of the TTS show variability in X-rays with 
time scales of minutes to days and can present strong flares. 

X-ray spectroscopy of CTTS has shown some evidence that 
their X-ray emission could be due, at least in part, to accre- 
tion shocks on the stellar photosphere (KastneretaT][2002a; 
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Stelzer & Schmitt 112004 ISchmitt et al. 112005b iGunther et al. I 
20061 lArgiroffi et al.l l2007h . although some CTTS do not 
support this me chanism. Their X-ray emission is believed 
to be coron al dAudard et al. ||200H iGudel & Telleschill2007l: 
LSmith et al. Il2005h . A X -ray s oft excess in CTTS is, h owever, 
reported compared to WTTS (Tell eschi et al. Il2007al) . X-rays 
in WTTS are thought to originate only from magnetic activity 
similar to main-sequence magnetically active stars, albeit at 
much higher levels. 

The low level of accretion in WTTS makes them the best- 
suited objects to study the magnetic activity in PMS stars, and 
the impact of stellar flar es and X-rays onto their optically thin 
disks and their planets (|Lammer et al. 1 12003; La mmer et al. I 
2006; Smi th & Scaloll2007l) . 

WTTS show s trong X-ray and non- thermal ra dio emissions 
dSkinner| [T993; Fei gelson & M ontmerle"l [l999[) . Compared 
to zero-age main sequence stars, their X-ray luminosities 
are very high (10 28 5 - 10 31 erg s" 1 ), but they typically show 
Lx/L bpi ~ 10~ 4 . WTTS d i splay very high cor onal tempera- 
tures (Skin ner et ainil997t ITsuboi et al. Ill998h . X-ray flares 
with temperatures of a few tens of MK a re freq uently seen, 
e.g., ITsuboi et al. I d!998l) ; IStelzer et al. I d2000l) . Thanks to 
the Chandra X-ray Observatory (CXO) and the X-ray Multi- 
Mirror Mission (XMM - Newton), observations of PMS stars 
such as the WTTS HDE 245059 can provide a better under- 
standing of magnetic activity and its evolution, allowing to 
resolve individual spectral lines and to get improved diagnos- 
tics of the plasma properties. 

CTTS and WTTS show similar ages and are found in the 
same area of the Hertzsprung-Russell diagram; why are they 
different then? The environment might be an interesting cause 
of the different properties. For example, a supernova explo- 
sion can clear away material, leading to the loss of the disks 
in some young stars (Dolan & Mathieu 1999, 2001, 2002). 
The presence of a massive star can also disperse the young 
circumstellar disks by photoevaporation due to UV radiation. 
On the other hand, the above cannot hold for low-mass star 
forming regions, such as Taurus. However, understanding the 
effect of environment onto young stars could give clues about 
the mechanism leading to the differentiation between CTTS 
and WTTS. 

The main goals that we address in this work are: i) to obtain 
the plasma properties; temperatures, abundances, X-ray lumi- 
nosity of the WTTS HDE 245059 located in the A Ori star 
forming region, ii) to obtain fluxes from individual lines, Hi) 
to derive age and masses of the system based on NIRSPEC 
data. 

2. THE HDE 245059 SYSTEM 

The young binary HDE 245059 is a member of the A Ori 
cluster, a star forming re gion located at d = 400 ± 40 pc 
dMurdin & Penston I fl977h . The A Ori cluster or Collinder 
69 includes the O star A 1 Orio nis, spectral type 08 III 
dBarrado v Navascue s et al.1l2007l) . The star is near the cen- 
ter of the region which is believed to be the remnant of a su- 
pernova that exploded about 1-2 Myr ago. This explosion 
might have cleared out the regi on, leaving a ring of mo l ecular 
gas s urrounding A 1 Orionis dDolan & Mathieulll999L 1200 ll 
2002). Star formation is, however, still active at the edges of 
the molecular ring. 

The A Ori region is also populated by faint PMS 
stars and has been studied a t different wavelengths; ra- 
dio, infrared, visible, UV dMaddalena & Morris I 1 19871; 
iDolan & Mathieu 1 119991; iBarrado v Navascues et a~f 120071) . 



Maddalen a & Morris! d!987l) argued that the observed CO 
emission from transition 7=1^0 together with UV, opti- 
cal, IR, and 21 cm data of the surrounding molecular clouds 
indicates that the ring centered on A 1 Orionis is actually the 
remnant of a preexisting cloud. Dolan & Mathieu (1999, 201, 
2002) have studied the region and attempted to deduce its his- 
tory. According to the authors, about 6 Myr ago the star for- 
mation process started in the most massive clouds of the re- 
gion giving birth to several OB stars. This process increased 
gradually until approximately 1 Myr ago when a supernova 
exploded disrupting the central region and decreasing the star 
formation rate, yet not stopping it completely. 

Our target (A Ori X- 1 ) was discovered in the X-rays with the 
Einstein Observatory during observations of the region cen- 
tered at A 1 Orionis. The photometric and astrometric study 
of the region shows evidence that our X-ray source is in fact 
member of the association. HDE 245059 shows a strong 
absorption in the Li (6707 A) line, an ind icator of youth 
dAlcala et alJ l2000h . iFernandez et al. I d 19951) has measured 
the equivalent width of the Ha emission line for HDE 245059, 
W(Ha) = 0.37 A, indicating that the star is a WTTS. The 
star has a spectral ty pe Kl, its X-ray lum inosity is high, 
logLx = 31.7 ergs s -1 dStone & Taam|[l985l) . Based on com- 
pariso n with evolutionary model isochrones, IStone & Taam I 
(1 19851) estimated a mass and age ofM = 2- 3M^ and? ~ 1-4 
Myr. More recent isochrones, e.g those bv lSiess et alJ d2000h 
might, ho wever, change t he abo ve estimates (see below, Sec- 
tion [5TTJ. Skinn er et al. I d!991l) obtained an upper limit esti- 
mate for the disk mass around HDE 245059 of 0.32 based 
on the 1 . 1 mm continuum emission. iPadgett I d 1996b obtained 
the effective temperature of HDE 245059, T eff = 54 1 ± 1 1 K, 
and an estimate of the photospheric iron abundance, [Fe/H]= 
-0.07 ±0.13. Our target is a fast rotator (v sin i ~ 25 km s -1 ), 
which in general is a sign of ste llar youth or bina ries of spec- 
tral types G and K dFekelll997|). The studies of lFekell d 19971) 
in the optical, and lAlcala et alJ d2000t) in the X-ray did not find 
evidence of radial velocity variability, an indication that HDE 
245059 is not a spectroscopic binary. However, our high spa- 
tial resolution near-infrared observations with Keck in dicat e 
that it is, in fact, a spatially separated binary (see Section lXTl ). 

3. OBSERVATIONS 

3.1. Near-infrared 

On Dece mber 13, 2003 we o bserved HDE 245059 with 
NIRSPEC dMcLean et al. 1 120001) installed behind the adap- 
tive optics system dWizinowic h et al. Il2000h on the 10m Keck 
II telescope. The imaging camera within NIRSPEC offers a 
pixel scale of 0"0168±0"0001 and an absolute orientation of 
l.°l±0.°8 as determined from observations of calibration bi- 
naries and of a reference field in the Orion Trapezium. HDE 
245059 was used as the adaptive optics guide star with the 
wavefront sensor running at a frequency of 488Hz, resulting 
in good correction (Strehl ratios of 43 and 19% at 2.2 and 1.6 
/im, respectively). We acquired images of the system with the 
broadband K (A = 2.20 /mi, AA = 0.39 /mi) and NIRSPEC-5 
(Ao = 1.61 /mi, AA = 0.40 /mi, a close analog to the usual 
H filter) filters as well as with the narrow band Br7 filter 
(Ao = 2.165 /mi, AA = 0.02 /mi). In each filter, a series of 
short exposures were coadded at four successive positions on 
the detector. The total integration times were 40s, 4s and 20s 
with the K, H and B1-7 filters, respectively. For each filter, 
the four independent images were medianed to create a sky 
that was subtracted from each image. The images were then 
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FIG. 1 . — Near-infrared Keck adaptive optics images of HDE 245059 in H, 
K, and Br7 (2.165 (im) bands. The image is in square root scale, North is 
toward the top and East is toward the left. 

corrected for flat-field effects and cosmetically cleaned. The 
images were then realigned and averaged to produce the final 
images. Relative astrometry and photometry was obtained us- 
ing the DAOPHOT package. The astrometric results from all 
three filters were finally averaged to reduce random uncertain- 
ties. 

The near-infrared images resolved HDE 245059 into a bi- 
nary. Figure Q] shows the near-infrared Keck adaptive optics 
images of HDE 245059 with the three filters. The northern 
component is brighter in all images. The projected separa- 
tion and uncertainty is 0."866 ± 0."005 and the position angle, 
measured East from North is 150.°0±1.°0. The final astro- 
metric uncertainties are dominated by the absolute calibration 
uncertainties of the detector. We measured magnitude differ- 
ences of 0.98 mag, 0.88 mag and 0.86 mag with the H, K, 
and Br7 filters respectively (with typical uncertainties of 0.03 
mag). These differences suggest that both components have 
similar broadband colors. 

3.2. X-rays 

X-ray observations were performed using the High-Energy 
Transmission Gratings (HETG) in combination with the Ad- 
vanced CCD Imaging Spectrometer (ACIS-S), on board of the 
Chandra X-Ray Observatory (CXO). We obtained a total ex- 
posure time of 93 ks scheduled in three epochs: December 
30, 2005, January 7, 2006 and January 13, 2006 (Observation 
identification numbers 6241, 7253, and 5420, respectively). 
The short time interval between the observations allowed sim- 
ilar roll angles: 308 .°5 for the first two and 294 .°8 for the last 
one. 

Both components of the binary were detected in the zeroth 
order image. The binary orientation was close to the disper- 
sion direction of the Chandra MEG arm. Choosing the origin 
of the wavelength system between the two stars we were able 
to separate them in the grating spectra despite the small sepa- 
ration. 

4. CHANDRA DATA REDUCTION AND ANALYSIS 

The data of the HETGS observations were reduced from 
level 1 event files with the Chandra Interaction Analysis 
of Observations software, CIAO 3.4, using the calibration 
database, CALDB 3.3.0.1, and following standard procedures 
to obtain a type 2 event file. We removed streak events which 
affect significantly ACIS-S4 (CCD ID=8). 

4.1. Zeroth Order Images 

For the zeroth order images we used the subpixel event 
repositi onary (SER) algorithm in order to improve spatial res- 
olution (iTsunemi et al. KOOltlEiet al. I2003ll2004h . When an 
X-ray photon hits the detector, the charge cloud created can 
either be spread into neighbouring pixels (called split pixel 
event) or remains in a single pixel (called single pixel event). 
Positions are determined with higher accuracy in the case of a 
split pixel event, however these events represent only a small 
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FIG . 2. — Zeroth order X-ray image of HDE 245059 in square root scale for 
the three observation epochs. A flare is visible in the south component during 
the first epoch. At the far right, the average image obtained by summing the 
threee observations. 

fraction of the total events. The SER technique uses both sin- 
gle pixel events and 2 pixel split events to increase the statis- 
tics. When applied to the ACIS observations the SER tech- 
nique reduces the uncertainties in the determination of the 
photon impact position improving the spatial resolution. The 
improvement in FWHM is typically between 40 to 70 %. 

The zeroth order images of HDE 245059 for each obser- 
vation epoch are displayed in Figure [2] A flare from the 
south component is visible in the first observation epoch (first 
panel). The summed zeroth order image over all the obser- 
vation runs, Figure [2] last panel, clearly shows the two com- 
ponents, the north component being on average the brighter 
one. 

4.2. Light Curves 

Figure [3] shows the light curves of the zeroth order data for 
the 3 observation epochs. The light curves were extracted 
from the data treated with the SER algorithm using a 2."1 ra- 
dius for the binary and 0"4 radius for each component, thus, 
the sum of the lightcurves of both components do not match 
the binary lightcurve. During the first observation epoch, a 
flare is detected from the south star, starting approximately 
6 ks after the beginning of the observation. The flare lasted 
for nearly 6 ks, peaking 1.5 ks after its start with 0.06 cts s , 
roughly 5 times higher than the southern star average count 
rate for the first observation epoch. 

We have also included the hardness ratio for the first ob- 
servation epoch as an inserted figure in the lightcurve plot 
(Figure [3] top panel). The hardness ratio was calculated us- 
ing the expression H/S, where the soft band (S) was taken in 
the range from 0.3 to 2 keV and the hard band (H) was taken 
in the range from 2 to 10 keV. 

4.3. Spectra 
4.3.1. Binary Components 

Thanks to the spatial resolution of Chandra we have sepa- 
rated the components of the HDE 245059 binary in the zeroth 
order spectrum. The spectrum was extracted after merging the 
event files from the three observation epochs following stan- 
dard CIAO procedures, using a circular region of radius 0"4 
for each star. The background was extracted from an annular 
region with radii of 4 "6 and 21" We have considered events 
within the range 0.2-7 keV, where most of the signal is con- 
centrated. We have obtained the spectra for the quiescent state 
for each component. 

In Figure [4] we have plotted the spectra and overlaid the 
best-fit model for each star. The northern star spectrum shows 
a higher signal with respect to the southern companion, which 
is consistent with the results from the zeroth order image. 

We have fitted the individual spectra of the binary compo- 
nents obtained from the zeroth order data (see Table [TJ us- 
ing XSPEC version 1 1 (lArnaudllT9 96). and a discrete emis- 
sion measure distribution (EMD). The hydrogen column den- 
sity and abundances were fixed to the zeroth order plus grat- 
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FIG. 3. — Zeroth order lightcurve of HDE 245059 for the three observation 
epochs, binned by 10 3 s. A flare is visible in the south component during 
the first epoch. The combined light curve of the binary is shown by the solid 
line, by the dotted line the contribution from the northern component and by 
the dashed-dotted line the contribution from the southern component. For the 
first epoch lightcurve we have also included a plot of the hardness ratio for 
the southern star. 



ing best-fit values (discussed in section |4. 3. 31 >. For the north 
component we have fitted a 3-T plasma, with temperatures 
in the range between 6 and 40 MK. The emission is domi- 
nated by the softer plasma with temperatures in the range be- 
tween 6 and 1 3 MK. The average temperature was defined as 
logr av = (S; log 7- x EMi)/EM tota i, leading to r av =11.3 MK. 
The upper limit for the component at 40 MK was not well 
constrained. 

For the southern star we could not fit the soft component 
at 6 MK that was obtained for northern star. We have found 
a 2-T plasma with temperatures 8 and 34 MK, the emission 
is dominated by the plasma at 8 MK. Again in this case, the 
upper limit for the highest temperature (34 MK) was not well 
constrained. The average temperature was 12.2 MK, slightly 
higher than the average of the northern star. 

We remark that the sum of the emission measures, as well 
as the luminosity, obtained for both components of the binary 
does not match the value of the total emission measure found 
in our fits to the combined spectrum of the binary (see !4.3.3"T l 
due to the different extraction radii used for the binary and 
their components. 



FIG. 4. — Zeroth order average spectrum over the three observation epochs 
(in energy) for each component of the binary HDE 245059 (north and south). 
We have overlaid the best-fit model (see Tablefl} as a solid histogram. 

4.3.2. Flare 

To analyze the spectrum from the southern star during the 
flare we have extracted only the zeroth order spectrum, since 
we do not have enough signal in the gratings. We have fitted 
a multi-T model leaving the photoelectric absorption compo- 
nent, abundances and temperatures fixed to the values from 
the best-fit to the quiescent state of the southern star. We 
have then added an extra temperature component for which 
we have left only the emission measure and temperature free 
to vary. We could not add more than one temperature compo- 
nent to our model, probably due to the low signal obtained. 

Our fit shows that the flare emission is dominated by a 
plasma at ~ 30 MK (see Table [2]). The luminosity during 
the flare increases to 7.6 x 10 30 erg s" 1 , which is a factor of 
five higher than the luminosity of the quiescent state. The 
high temperature found during the flare is also present in the 
quiescent state, indicating that the variable component during 
the flare originates mainly in the hotter part of the emission 
measure distribution. 

4.3.3. HDE 245059 Binary 

The combined zeroth order CCD spectrum of HDE 245059 
was extracted after merging the event files from the three ob- 
servation epochs following standard CIAO procedures. This 
was possible thanks to the similar pointing and setup between 
the observations. We have used a circular region with 3 "9 ra- 
dius. The background was extracted from an annular region 
with radii of 4"6 and 21" We considered events at energies 
within the range 0.2-7 keV. 

The grating spectra of the HDE 245059 binary were ex- 
tracted for both the HEG (High Energy Grating, 0.8-10 keV) 
and MEG (Medium Energy Grating, 0.4 - 5 keV) orders ±1. 
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We have selected the wavelength range to 2.8 - 25 and 1.8 - 
17.3 A for MEG and HEG, respectively. The spectra were 
binned by a factor 3 in wavelength, in particular to sum the 
contribution of both stars in the MEG spectra (since both 
stars are aligned with the dispersion direction). The respective 
spectral responses were generated using standard CIAO tools. 
The spectra of all the observation epochs were then merged 
using the CIAO procedure merge_all. In Figure [5] we show 
the binary grating spectrum for MEG and HEG first orders. 
We have labeled some important emission lines. 

We have analyzed the spectrum of the binary HDE 245059 
during the quiescent state using the combination of the zeroth 
order spectrum plus the grating spectrum from HEG order 1 , 
and MEG order 1. For the fits we have used three methods: 
a multi-T component model as a discretization of the EMD, 
a continuous EMD form Chebychev polynomials, and a con- 
tinuous EMD approximated by two power laws. For the three 
methods we have obtained abundances, emission measures, 
and an estimate of the photoelectric absorption component 
parameterized by the hydrogen column density. In all cases 
we have used the combination of the summed zeroth order, 
HEG±1 and MEG±1 spectra. 

4.3.4. Method 1: Discrete Emission Measure Distribution with a 

Multi-temperature Model 

This method uses a discretization of the EMD by a multi- 
temperature VAPEC model, which calculates both line and 
continuum emissivities for a hot, optically thin plasma. We 
have combined this model with a photoelectric absorption 
component (WABS), parametrized by the neutral hydrogen 
column density Nr. R eference abundances were set to so- 
lar photospheric values dGrevesse & SauvaT1ll998l) . We have 
used several isothermal components for the plasma emission. 
A 4-temperature plasma model gave the best fit compared to 
models with additional or fewer plasma components. Our re- 
sults are summarized in Table [3] All the values are displayed 
with their respective errors calculated for a confidence range 
of 68%. We have found a plasma between 4 and 50 MK. 
The average temperature was defined as logr av = (£,log7; x 
EMO/EMtotai, leading to T av = 10.7 MK. The value of the hy- 
drogen column is very low, Nu = 8 x 10 19 cm -2 , corresponding 
to Ay = 0.04 mag, which is consistent with the 2MASS col- 
ors of the binary (J-K = 0.6) and with photospheric colors of 
early K type stars. We remark that the results obtained from 
the fit to the quiescent state are not different from the results 
obtained from fitting the total average spectrum (considering 
also the flare). 

The best-fit model overlaid to the HEG±1 and MEG±1 
data for the quiescent state are shown in Figure[5] The bright- 
est line detected is Ne X Lya at 12. 13 A. Significant emission 
comes also from O VIII Lya at 18.97 A, Ne X Ly/3 at 10.23 A. 
Interestingly there is no detection of O VII or N VII despite 
the low Nu- Lines from highly ionized states of Fe XVII to 
Fe XXIV confirm the presence of a wide range of plasma tem- 
peratures revealed by our fits. 

4.3.5. Method 2: Continuous EMD from Chebyshev Polynomials 
The differential EMD tp{T) is given by 

<p(T) = n H n e ^- (cm^K" 1 ) (1) 
dT 

where the total EM is given by EM tot = / ip(T)T AlogT 
= J (p(T)Td(lnT). A graphical representation of the EMD 



independent of the grid bin size (AlogT) is given by 
EMD(r)/Alog7\ This method assumes that the shape of (p(T) 
can be approximated by the exponential of a polynomial given 
by ip(T) = ae^ T \ where a is a normalization constant and 
ui(T) is a polynomial function of the temperature, which we 
have chosen to be a Chebyshev polynomial (see Lemen et al. 
1989, Audardetal. 2004). 

Our model uses a grid of temperatures in the range 
logT(K) = 8-10 with AlogT = 0.2 dex for a polynomial de- 
gree of n = 8 which gives the optimal fit. Coronal abundances 
and the photoelectric absorption were left free to vary. The 
results obtained with this method are displayed in Table|4] 

4.3.6. Method 3: Continuous EMD approximated by two power 

laws 

This method is based on a continuous EMD described by 
two power laws; one at low temperatures with slope a, and 
one at high temperatures with slope (3 (see Telleschi et al. 
2007b). The EMD peaks at the temperature To which also 
represents the limit between the two regimes. 

arn - / EM o(T/To) a for T < T 
^ ' {EM (T/T f for T>T 

The normalization parameter is defined as the EM in the 
temperature bin at 7J). The free parameters are: 7b, EMq, a, (3, 
Nu, and the elemental abundances. We have used 2 approx- 
imations; first leaving a and f3 free to vary, and after fixing 
the value of a = 2 leaving (3 free to vary. The value a = 2 is 
usually found in magnetically active main-sequence and pre- 
main sequence stars. Alogr=0.1 dex was set in both cases. 
The results are presented in Table [5] 

In Figure|6]we have plotted the EMD of the binary obtained 
by the different methods applied. In order to compare them 
we have used the quantity EMD(r)/Alogr which is indepen- 
dent of the bin size. The total volume emission measure ob- 
tained from the different methods are similar: 7.3 for method 
1, 7.1 for method 2, 7.5 for method 3 a fixed and (3 free, and 
7.3 for method 3 with a and j3 free (all values are in units of 
10 54 cm" 3 ). 

4.4. Ftting 
AAA. Electron Densities 

We have obtained the fluxes of the density sensitive He-like 
triplets Si XIII, Mg XI, and Ne IX from the binary grating 
spectrum (HEG and MEG first order). The first step was to 
fit the continuum to get an estimate of the plasma tempera- 
ture. We used a bremsstrahlung model, taking the spectrum 
without the contribution of the bright lines and considering 
the wavelength range between 1.2 and 40 A. We have ob- 
tained a plasma temperature of 12.4 ±0.3, MK close to the 
average temperature found in our previous best-fit of the grat- 
ing plus zeroth order spectrum. To fit the triplet lines (res- 
onance (r), intercombination (/), and forbidden (/)) we used 
the combination of delta functions for the line profiles and the 
above mentioned bremsstrahlung model for the continuum. 
We have fixed the continuum parameters and fitted the triplets 
only in the wavelength range of interest. The only free param- 
eters in our fits were the line fluxes, since the wavelength of 
each line was fixed. Since the Ne IX triplet is blended with 
Fe XIX, we added a delta profile to consider this blend. The 
results are summarized in Table [6] In Figure [7] we have plot- 
ted the three triplets for the MEG first order spectrum includ- 
ing the delt a profiles used to fi t the lines. Based on theoreti- 
cal models dPorquet et al. Il200lb . using the R = f/i ratio and 
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FIG. 5. — HEG ±1 and MEG ±1 average spectra and best-fit model overlaid for the binary. Bright lines are labeled. The spectra were binned by a factor of 3. 





log T (K) 



FIG. 6. — Emission measure distribution of the binary HDE 245059 spectra 
obtained from the different methods described in the text. We have plotted the 
emission measure per bin; EMD(Y)/A log T, which is independent of the bin 
size allowing comparison with the other methods. In the case of the multi-T 
approach we have plotted the EM for each component(filled circles). For the 
upper panel we have used a linear scale in order to emphasize the dominant 
plasma at 10 MK. In the lower panel we have used logarithmic scale to show 
the presence of weaker components of plasma at very high temperatures. 



the temperature of the emitting plasma we have obtained an 
estimate of the plasma electron density. We have calculated 
the densities using confidence levels of 68 and 90% (see Ta- 
ble |7J. For the 68 % confidence level, we have obtained n e = 
5 ±5 x 10 1 1 cm" 3 from the Ne IX triplet, n e = 5 x 10 13 cm" 3 
from Mg XI, and n e < 5 x 10 13 cm" 3 from Si XIII, which is 
consistent with the low density plasma. For the 90 % con- 
fidence level, we have obtained n e =< 2.0 x 10 12 cm" 3 from 
the Ne IX triplet, n e = 1^ 8 x 10° cm" 3 from Mg XI, and 
n e < 3 X 10 14 cm" 3 from Si XIII. Therefore, we prefer to err 
on the safe side and conclude that there is no evidence of high 
densities in HDE 245059. 

4.4.2. Lines fluxes from each binary component 

We have also attempted to obtain individual fluxes from 
several lines for each component of the binary. We have used 
only the first order spectra of MEG because of the higher 
signal to noise, and because the binary was aligned along 
the MEG dispertion direction, allowing us to disentangle the 
two components in wavelength space. The spectrum was re- 
extracted using a bin size of 0.015 A. To fit the single lines 
we used the method described in section l4.4.11 a combination 
of a delta function for the line profile and the bremsstrahlung 
model for the continuum. The fits to the MEG+1 and MEG-1 
data were made simultaneously because the shift of the lines 
will occur in opposite directions from the origin of the ref- 
erence wavelength. For each line we considered only the in- 
strumental profile, therefore we had 4 delta profiles: 2 for 
the MEG+1, northern and southern star and the same for the 
MEG-1. From the grating equation, we calculated the shift 
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of the lines. The line shifts were smaller than the line spread 
function full width at half maximum (FWHM), with a typi- 
cal separation of 17 mA. To fit the lines, we fixed the energy 
at which they were expected to be found, leaving only the 
fluxes free to vary. Besides, the fluxes from each star were 
linked for MEG+1 and MEG-1. The fits thus considered 2 
free parameters: the flux of the line from each star. The cal- 
culations were made whenever the signal for a given element 
was high enough; we have then obtained fluxes for O VIII 
Lya, Ne X Lya, Ne X Ly/3, and Mg XII Lya. For the He-like 
triplets the signal was too low, thus we did not obtain fluxes 
for the single binary components. The results are displayed in 
Table [8] in Figure [8] we have plotted the line profiles for the 
full resolution grating MEG± 1 data. 

5. DISCUSION 

5.1. Environment and Binarity 

The absence of a disk in the WTTS HDE 245059 
members despite their young age might be closely re- 
lated to the environment where this young binary has been 
formed. The evolutio n of the region around A Ori has 
been discussed in detail dDolan & Mathieu 1 1999] 1200 ll 120021: 
iBarrado y Navascues et al. 1120071) . Dolan & Mathieu (1999, 



Si XIII 




6.5 6.6 6.7 _ 6.8 

Wavelength (A) 



MgXI 




9.1 9.2 9.3 

Wavelength (A) 




13.3 13.4 13.5 13.6 13.7 

Wavelength (A) 



FIG. 7. — He-like triplets in the grating spectra. We have plotted the MEG 
spectrum of the binary and in histogram the delta profiles used to fit the res- 
onance (r), intercombination (i), and forbidden (/') lines. For Ne IX we have 
also added the blend with Fe XIX. 



2001, 2002) presented the hypothesis of a supernova explo- 
sion that cleared out the region about 1 Myr ago leaving a 
molecular ring. During the phase prior to the supernova ex- 
plosion the stars in the cluster would have been confined in 
a small region and the circumstellar disks cleared away by 
photo-evaporation. This hypothesis is sup ported by the small 
fracti on of CTTS in the region, ~ 7 % dDolan & Mathieu I 
1999), whic h is low when compared with cluste rs with similar 
properties. IBarrado v Navascues et al. 1 (120071) have studied 
the members of the A Ori with the Spitzer satellite, finding a 
31% of members with disks in A Ori, but only 14% with thick 
disks. According to this study, the presence of CTTS near the 
center of the cluster would suggest that massive stars and the 
supernova explosion had no major effect on the disks, or that 
they have been formed after the explosion. Rece nt high reso- 
lution optical spectroscopy of the A Ori cluster (Sacc o et al. I 
2008) has given a fraction of stars with disks of 28%, higher 
than the previous studies. Our new estimates of the properties 
of the HDE 245059 binary (see section 15.2b give an age of 
w 2-3 Myr: at this age the binary would have experienced the 
supernova explosion, according to Dolan & Mathieu (1999, 
2001, 2002) hypothesis. Another aspect to be considered is 
the influence of binarity i n the disk truncation . This hypothe- 
sis has been discussed by Kra us et ai~l ([2008), who presented 
recent results showing that several of the young stars with- 
out disks in a survey of nearby star-forming regions are close 
binaries. 

5.2. Stellar properties 

The discovery of the binarity of HD 245059 prompts us to 
re-evaluate the stellar properties of the system. The few phys- 
ical properties; mass, age, spectral type, and effective tem- 
perature we re previously obtained under the assu mption of a 
single star (Pad gett Ill996t IStone & Taamlll985l) . Our near- 
infrared and X-ray data have been the first to separate it into a 
binary. The question that arises immediately is how different 
the properties of the HDE 245059 members are. The mag- 
nitude differences from the near-infrared images in the H, K, 
and Br7 bands suggest that both components have similar col- 
ors, but the northern star is brighter in both near-infrared and 
X-ray images, which might indicate that is the more massive 
of the system. Throughout this analysis, we assume a distance 
of 400 pc to the system. We further assume that the flux we 
receive from each component can be entirely attributed to the 
stellar photosphere, without contribution from accretion for 
instance, in line with the WTTS status of the system. 

In a first approach, the absolute near-infrared magnitudes 
of each component can be used to estimate its main proper- 
ties (mass and age). Combining the unresolved 2MASS pho- 
tometry of the system with our new H and K band flux ra- 
tios, we determine absolute magnitudes of Mh = -0.07 and 
M K = -0.20 for the primary and M H = +0.91 and M K = +0.68 
for the secondary 10 . Such absolute magnitudes are too bright 
for solar-like stars even at ages as young as 1 Myr, base d on 
the stellar evolut ionary models of Bar affe et al. I (119981) and 
ISiess et ail d2000l). To exp l ore hig her-mass regimes, we adopt 
the models of ISiess et al.l (120001) which extend up to 7M . 
Based on this model, for ages of 4 Myr or more, only B-type 
stars reach the observed brightness of HD 245059 north and 
south. At 3 Myr, at least the primary would have to be a B 
star. Since this can be confidently excluded from the spectral 

10 Note that extinction, if present in front of the system, would yield even 
higher absolute magnitudes. 
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FIG. 8.— The plot shows the Mg XII Lya, Ne X Lya, Ne X Ly/3, and O ' 
MEG+1. The line were modeled with a bremsstrahlung plus a gaussian mode 
purely instrumental. 

analysis of Padgett ( 1996), we conclude that the system is no 
older than 2 Myr. Assuming stellar ages in the 1-2 Myr range 
(since st ar formation appea r s to ha ve ceased about 1 Myr ago, 
see e.g. iDolan & Mathieu I ( 1200 lb ), a primary of 2.5-3.5M Q 
and a secondary of 2-3 M Q would match the observed near- 
infrared magnitudes of the two components. 

To improve on this estimate and take advantage of a broader 
dataset, we perform a fit to the optical and infrared SED of the 
system. For this purpose, we use t he unresolved UBV pho- 
tometry from lStone & Taarnl (Tl985). the unresolved J magni- 
tude from 2MASS, the spatially resolved photometry derived 
above and the 3.6, 4.8, 5.6 and 8 fim unresolved IRAC pho- 
tometry determined from archival images and using default 
recipes for aperture correction around point sources. This 
represent a total of 12 independent me asurements. We used 
a grid of NextGen stellar spectra from Bar affe et al. I (Il998h 
with logg = 4.0 (our results are largely insensitive to the as- 
sumed surface gravity). We used 5 free parameters in our 
model: the stellar effective temperatures T e ff(N) and T & ff(S), 
the stellar radii and R$, and the extinction Ay. To en- 
sure that we did not miss the best possible solution, we con- 
servatively allowed for wide ranges of initial guesses (4600- 
7600 K for r eff (A0, 2400-7000 K for r eff (S), 3.5-6 R Q for R N , 
2-9 Rq for Rs and 0-1.5 mag for Ay). In this procedure, we 
do not use information from evolutionary models as prior and 
only require that T e ff(S) < T e ff(N). Overall, we tested about 
5 million independent combinations of the 5 free parameters, 
estimating a reduced \ 2 value for each model. We then used 
a Bayesian inference method to explore the parameter space: 



Til Lya lines for the grating data. In solid line the MEG- 1 and in dashed line the 
taking into account the contribution of both stars. We considered the profiles to be 

each model in the grid is assigned a probability p = e~ x ~l 2 , and 
1- and 2-dimensional probability distributions for individual 
and pairs of free parameters are then produced by marginaliz- 
ing the hypercube against the other dimension. 

Using the mode of each 1 -dimensional probability distri- 
bution and defining a 68% confidence level interval around 
it, we infer T eff (N) = 5850^ K, T eH (S) = 3460^° K, R N = 
4.94^-^Rq and R s = 4.27+^/?©. However, as illustrated in 
Figure |9] there is a substantial ambiguity between T e ff(S) and 
Rs due to the fact that we only have 2 measurements that di- 
rectly constrain the secondary. A cool but large secondary 
fits equally well the data than a warmer and more compact 
star. This correlation between the stellar parameters also ex- 
plains the difference between the absolute best model (i.e., 
lowest x 2 ) and the most probable stellar parameters based on 
the 1-dimensional probability distributions. Both estimates 
nonetheless agree within the uncertainty, although we prefer 
to use the 2-dimensional probability distributions to estimate 
the stellar properties. 

The final step in this analysis co nsists in overplotting the 
predictions of the lSiess et al.l d2000) evolutionary models. We 
note that both the stellar luminosity and radius are direct out- 
put of these models so this is a natural set of parameters to 
compare model and data. We overplot in Figure |9]the 1, 2, 3 
and 5 Myr isochrones and readily conclude that the primary 
star is most likely w 2 Myr-old. Similar to the conclusion 
based solely on the near-infrared magnitudes, stellar ages be- 
yond 3 Myr can be confidently excluded. While the best fit 
to the secondary suggests an age even younger than 1 Myr, 
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the 2 Myr isochrone does intercept the 68 % confidence level 
contour. Using 6000 K and 5Rq (equivalenfly , 29 L©) for the 
primary, interpolation in the lSiess e t al. (2000) model yields a 
stellar mass of 2.1 Mq and an age of 2.7 Myr. From the extent 
of the 68% confidence level contour, we estimate uncertain- 
ties of 0.5 Mq and 1 Myr. Based on the 2.7 Myr best-fitting 
age, we further infer M s = 2.3-2.4M© . 

While relatively large uncertainties remain due to the lim- 
ited number of resolved photometric measurements of the 
system, we conclude that the stellar masses are w 3M© and 
« 2.5M© and the age of the system is « 3 Myr. High spa- 
tial resolution optical data, which are currently unavailable 
for the system, would greatly improve the accuracy on these 
parameters. Nonetheless, we consider that this is a significant 
improvement over the previous estimates that did not take into 
account the fact that the system is indeed a binary. As a final 
note regarding the stellar properties, we note that the age of 
HD 245059 is younger than the average age of other members 
of the AOri association but also consistent with that of the 
youngest systems in the associations. 



8 -\\ 
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FIG. 9. — Estimated stellar radii and effective temperatures based on our 
fit to the SED. Solid and dotted contours represent the 68.3, 95.5 and 99.7% 
confidence levels for the primary and secondary, respectively. Filled dia- 
monds indicate the best fitting stellar models in our grid (corresponding to 
Ay = 0.7 mag). Open diamonds represent the most probable models from 
the mode of the 1 -dimensional probability distributions in our Bayesian ap- 
proach; associated uncertainties represent the 68.3% confidence level inter- 
vals around these models. From top to bottom, the dashed line s represent 
the 1, 2, 3 and 5 Myr isochrones from the evolutionary model of jSiess et aT] 
2000). The available dataset is consistent with a 2-3 Myr age for the system 
(see text for more detail). 



5.3. Spectral X-ray Properties 

The fits to the high-resolution grating spectroscopy data 
have revealed that both stars have similar spectral properties 
and emission measure distributions. Indeed, the spectra of the 
single components are consistent with the average spectrum 
of the binary. The similar temperatures of both stars allowed 
us to fit easily the grating spectra together. 

We do not detect the oxygen triplet in HDE 245059, ei- 
ther in the average or in the single spectra, despite the low 
TVh value. This triplet is used as an electron density indicator 
in the cool plasma component, expected to be present in the 
spectra of accreting stars. On the other hand we do find a soft 
plasma component at 3.8 MK in the average spectrum: the 
He-like triplet of Ne IX was detected with low signal to noise 



and it was found to be blended with the Fe XIX line. We also 
detect emission from Fe XVII blended with other lines. 

We have estimated an upper limit flux for the O VII triplet 
(see Table |6]l at 2 confidence levels: 68 % and 90 % for the 
average binary spectrum. 

In order to determine the nature of our non-detection of 
the O VII triplet, we have fitted some single lines for the 
combined spectrum of the HDE 245059 binary: O VIII Lya, 
Ne IX, Ne X Lya, and Fe XVII at 15 A (Table©. We have 
compared the results with the obse rved line fluxes of 3 young 
active st ars: 47 Cas B, EK Dra dTelleschi et al. 1120051) and 
AB Dor dGarcfa- Alvarez et al. Il2008h . We have found the line 
fluxes of the HDE 245059 binary to be about a factor of 100 
larger than the fluxes of the comparison stars. Using this ra- 
tio for the O VII resonance line at 90% confidence level, our 
upper limit for HDE 245059, < 2 x 10" 14 (ergs cm -2 s _1 ), re- 
mains close to what could have been expected based on the 
comparison stars. We conclude that the non-detection of the 
O VII triplet is probably due to a lack of sensitivity around 
22 A. 

Observations of the Sun show that abundances are related 
to the first ionization potential (FIP) of the elements in such a 
way that elements with low FIP (< 10 eV) are overabundant 
in the solar corona when compared to the photosphere and 
high FIP elements (> 10 eV) have similar ab undances in the 
corona and the photosphere (Feldman 1992). Observations 
of young, magnetically active stars show that the FIP effect 
is inversed with respect to the Sun, i.e. low FIP elements are 
underabun dant in the stellar coronae relative to elemen ts with 
high FIP (iBrinkman et al. Il200lt lAudard et al. 1 12003b . Fig- 
ure[lO]shows the coronal abundances with respect to the solar 
photospheric values for our binary plotted against the FIP. The 
abundances follow the trend of an inverse FIP effect, except 
for the Ca abundance which has a higher value, but we recall 
that the uncertainty in Ca abundance from our results is also 
high. 
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FIG. 10. — Coronal abundances with respect to the solar photospheric val- 
ues, obtained by fitting the combined zeroth order, HEG, and MEG data of 
the binary. The plot shows the best-fit values using 4 isothermal plasma com- 
ponents plus interstellar absorption. We see evidence of an inverse FIP effect 

A common problem found in the studies of the coronal 
element abundance is the lack of measurements of the stel- 
lar photospheric abundances, as they are difficult to obtain 
due to the large rotational velocity of magnetically active 
stars. Furthermore, a good knowledge of the stellar param- 
eters and atmospheric models are also needed. To sort out 
this problem, the stellar photospheric abundances are used 
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as a reference set. In the case of HDE 245059 Padgett 
(1996) obtained the photospheric iron abundance [Fe/H] = 
-0.07 ±0.13 usi ng as reference s et the solar photospheric 
abundances f rom iGrevessel d!984l). U sing the revised ref- 
erence set of iGrevesse & Sauvall (1 19981) . the stellar photo- 
spheric abundance is [Fe/H] = +0.10. Our coronal abun- 
dances have been obtained from the best fit to the zeroth 
order, MEG±1, and HEG± 1 using as reference set t he so- 
lar photospheric values from Grevesse & Sauvall d 1998b . We 
have obtained an iron abundance of [Fe/H] = -0.64. Thus, 
the difference between the photospheric and coronal values is 
[Fe/H] photosp heric- [Fe/H] coronal = +0.74, i.e., the coronal iron 
abundance is 5.5 lower than the photospheric one. This result 
is consistent with the inverse FIP effect expected in young 
magnetically active stars. 

6. SUMMARY & CONCLUSIONS 

We have obtained the Chandra high resolution spectrum for 
HDE 245059. Thanks to our X-ray and near-infrared data we 
have resolved this X-ray luminous WTTS to be a binary. We 
have attempted to get an estimate of the properties of the sin- 
gle components of the HDE 245059 binary based in the com- 
bination of the infrared magnitude differences, evolutionary 
models and SED determination. Our analysis gave a system 
of ss 2-3Myr with masses M N « 3M Q and M s « 2.5M©, 
rather high values for low-mass pre-main sequence stars. This 
is the first attempt to constrain the binary component masses, 
further studies are needed to obtain more accurate estimates. 

In the X-rays we were able to resolve both binary compo- 
nents in the zeroth order image and in the grating spectra. We 
have analyzed the zeroth order spectrum for each component 
of the binary using a multi-T plasma model. For the northern 
star we have found a temperature between ~ 6 to ~ 40 MK, 
the dominating component being a plasma between ~ 6 and 
~ 13 MK. For the southern star, we have found temperatures 
between ~ 8 and ~ 34 MK. 

During our observations, split in three runs, we have de- 
tected a flare from the southern star, the fainter of the system 
in average. An analysis of the zeroth order spectrum of this 
star during the flare has gave a plasma with a temperature of 
~ 30 MK, consistent with the high temperature component 
obtained for the quiescent state, and a luminosity higher by a 
factor of ~ 5 when compared with the quiescent state. 

We have derived the properties of the plasma from the com- 
bined zeroth order, MEG±1, and HEG±1 data for the binary 
during the quiescent state. We have analyzed the spectrum 
using 3 methods: a continuous emission measure distribution 
from Chebyshev polynomials, a continuous emission measure 
distribution approximated by two power laws, and the classi- 
cal discretization of the emission measure distribution by a 
multi-T optically thin plasma model. The three methods are 
consistent and we have found the emission to be dominated 
by a plasma in the temperature range 8-16 MK, with a low 
value for the column density, Nh = 8 x 10 19 cm" 2 . The low 
Nh might be due to the clearing of the inner region of A Ori 
due to a supernova event. The X-ray luminosity of the binary 
is high, Lx ~ 10 31 ergs s" 1 . 

The abundance pattern shows an inverse FIP effect, except 
for Ca which could not be well constrained by our fits. This 
result was confirmed for Fe by comparison of the coronal 
abundance from this work with the photospheric value from 



Padgett (1996). Despite the low N H , we did not detect the 
density sensitive O VII triplet, but we did obtain an estimate 
of the flux upper limit for this triplet from the binary grating 

TABLE 1 

BEST FIT TO THE ZEROTH ORDER SPECTRA OF THE HDE 245059 
BINARY COMPONENTS 



Parameter 


north 


south 


Ti (MK) 






T 2 (MK) 


12 7 +L8 

39 7+ 2 " 
■"•'-8.9 


8.3+0-6 

m 7+23.2 
■"•'-11.6 


T 3 (MK) 


T av (MK) 


11.3 


12.2 


EMi (10 54 cnr 3 ) 


1 o +0 - 2 




EM 2 (10 54 cnT 3 ) 


0.9 + 0.2 


0.6 + 0.01 


EM 3 (10 54 cnT 3 ) 


4+0.1 


0.2 + 0.01 


N H (10 I9 cirT 2 ) 


:=7.7 


:=7.7 


Flux (10~ I3 ergs crrr 2 s" 1 ) 


1.7 


0.7 


L x (10 30 ergs s" 1 ) 


3.3 


1.4 


C-statistics 


94 


56 


dof 


99 


51 



NOTE. — We have included an absorption model with the value of the col- 
umn density fixed to the best fit model to the quiescent state for the combined 
zeroth order plus grating spectra of the binary. 



spectrum. We have compared the fluxes obtained from our fits 
to O VIII, Ne IX, Ne X, and Fe XVII with observed fluxes of 
three active stars. Our upper limit flux for the O VII triplet is 
consistent with the flux differences when compared with the 
active stars. Based on the 90% confidence range, we have ob- 
tained upper limits to the plasma electron densities from He- 
like triplets; n e < 2 x 10 12 cm" 3 for Ne IX, < 6 x 10° cm" 3 
for Mg XI, and < 3 x 10 14 cm" 3 for Si XIII. 

According to our analysis, the properties of the 
HDE 245059 binary in the X-rays are similar to what is 
observed in other WTTS. Its accretion history might have 
been affected by the clearing of the region surrounding A Ori, 
probably by a supernova explosion, but its coronal properties 
have not been strongly modified. 

M. Audard and C. Baldovin Saavedra acknowledge the sup- 
port of the Swiss National Science Fundation grant PP002- 
1 10504. Support from the Chandra award SAO GO5-6012X 
is also acknowledged. We would like to thank Joel Kastner 
and David Huenemoerder for useful discussion on the SER 
method, and the anonymous referee for useful comments. 

Data presented herein were obtained at the W. M. Keck Ob- 
servatory, which is operated as a scientific partnership among 
the California Institute of Technology, the University of Cal- 
ifornia and the National Aeronautics and Space Administra- 
tion. The Observatory was made possible by the generous 
financial support of the W.M. Keck Foundation. The authors 
wish to recognize and acknowledge the very significant cul- 
tural role and reverence that the summit of Mauna Kea has 
always had within the indigenous Hawaiian community. We 
are most fortunate to have the opportunity to conduct obser- 
vations from this mountain. Part of this work has been sup- 
ported by the National Science Foundation Science and Tech- 
nology Center for Adaptive Optics, managed by the Univer- 
sity of California at Santa Cruz under cooperative agreement 
AST 98-76783 and by the Packard Foundation. 



REFERENCES 



Alcala, J.M., Covino, E., Torres, G., Sterzik, M. R, Pfeiffer, M. J. & Argiroffi, C, Maggio, A. & Peres, G. 2007 A&A, 465L, 5A 

Neuhauser, R. 2000, A&A, 353, 186 



HDE 245059: A binary revealed by Chandra & Keck 



11 



TABLE 2 

Best fit to the zeroth order spectrum of HDE 245059 south 

DURING THE FLARE. 



Parameter 


value 


T(MK) 


30 4 +u - b 


EM (10 54 cm 4 ) 


If® 


N H (10 19 cnr 2 ) 


:=1.1 


Flux (10 - 13 ergs cm 4 s~') 


3.9 


L x (10 30 ergs" 1 ) 


7.6 


C-statistics 


8 


dot" 


11 



NOTE. — The fit was obtained fixing all the parameters: temperatures, 
emission measures, abundances, and hydrogen column density to the values 
obtained from the best fit to the quiescent state and adding an extra tempera- 
ture component which was set free to vary. 



TABLE 3 

Best fit 4-T model of the quiescent spectrum of the 
HDE 245059 binary using Method 1 



Parameter 


Quiescent 


Ti (MK) 


3.9 ±0.3 


T 2 (MK) 
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T 3 (MK) 




T 4 (MK) 


50.2f, 


T av (MK) 
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EMi (10 54 cm" 3 ) 


1.1 ±0.2 


EM 2 (10 54 cm" 3 ) 


2.9 ±0.3 


EM 3 (10 54 cm" 3 ) 


9 7+0.4 
-0.3 


EM 4 (10 54 cm- 3 ) 


6 +0 2 


EM total (10 54 cm- 3 ) 


7.29 


Nu (10 19 cnr 2 ) 


7.7+2.1 


O 


30+ '^ 5 


Ne 


0.71 ±0.06 


Mg 


0.21 ±0.03 


Al 


0.34!° : 2 


Si 


0.18 ±0.02 


S 


1 4+005 


Ar 


0.46 ±0.3 


Ca 


0.58 ±0.3 


Fe 


0.23 ±0.02 


Ni 


:=Fe 


tx (10-' 3 ergs cm -2 s -1 ) 


5.2 


L X (10 31 ergs- 1 ) 


1.0 


C-statistics 


4148 


dof 


3708 



NOTE. — Method 1 use s a discretization of the EMD by a multi- 
temperature model (see |4.3.4"l for details.) We present here the best fit 4-T 
model of the average HDE 245059 binary for the quiescent state. We have 
used the combination of the zeroth order data plus the gratings HEG± 1 and 
MEG± 1 . Abundances are given with respect to the solar photospheric values 
iGrevesse & Sauval 1998). All errors are calculated with AC=1. 
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TABLE 4 

best fit to the quiescent spectrum of the hde 245059 binary 
using Method 2 



Parameter 


value 


EMtotaj (10 M cm 4 ) 


7.3 


Nyi (10 19 cm- 2 ) 


6.8 + | 


O 


0.34 ±0.03 


Ne 


0.69 ±0.06 


Mg 


0.21 ±0.02 


Al 


0.33 ±0.2 


Si 


0.18±0.02 


S 


1 5+ 005 


Ar 


0.46 ±0.3 


Ca 


0.51 ±0.3 


Fe 


0.23 ±0.02 


Ni 


= Fe 


C-statistics 


4161 


dof 


3708 



NOTE. — Method 2 is based on a continuous EMD obtained from Cheby - 
chev polynomials, in this case we have used a degree of 8 (see 14.3.5) . 
This method has been applied to the quiescent state spectrum of the bi- 
nary. Abundances are giv en with respect to the solar photospheric values 
IGreves se & Sauval 1998). Errors are calculated with AC=1. 



TABLE 5 

best fit to the quiescent spectrum of the hde 245059 binary 
using Method 3 



Parameter 


a fixed 


a, [3 free 


a 


:=2.0 


2.4 ±0.4 


P 


-1.5 ±0.09 


-1.5 ±0.09 


log To 


7.0±0.01 


6 n+0.03 
o - y -0.02 


EM tota i (10 54 cm- 3 ) 


7.49 


7.31 


Nh (10 19 cm- 2 ) 


8.7 ±2.5 


8.0 +2 « 


O 


0.28 ±0.04 


30+° 05 
-QQ4 


Ne 


0.68 „„ 

21 +i > 03 
u - -0.02 


69+™ 


Mg 


20 +i > 03 
u.zu 002 


Al 


0.35±0.2 


0.34 ±0.2 


Si 


0.18 ±0.02 


0.17 ±0.02 


S 


15+ 005 

u - 1J -0.04 


15+ 005 

u - 1J -0.04 


Ar 


0.43 ±0.3 


0.44 ±0.3 


Ca 


0.45 ±0.3 


0.46 ±0.3 


Fe 


0.23 ± 0.02 


0.22 ±0.02 


Ni 


= Fe 


= Fe 


C-statistics 


4168 


4167 


dof 


3714 


3713 



NOTE. — Method 3 is based on a continuous EMD described by two 
power laws. It has been applied to the quiescent spectrum of the binary. In 
the first case we have fixed one of them, whil e in the second case we have left 
both the power law indexes free to vary (see |4.3.6l for d etails). Abundances 
are gi ven with respect to the solar photospheric values iGrevesse & Sauval 
1998) Errors are calculated with AC=1. 
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TABLE 6 

Line fluxes for the binary HDE 245059 



Ion 


A (A) 


Flux (10 5 photons cm 2 s ') 


Si XIII(r) 


6.65 


0.5 ±0.1 


Si Xffl(i) 


6.68 


0.1 ±0.04 


Si XIII(f) 


6.74 


0.3 ±0.1 


Mg XI(r) 


9.17 


0.3±0.1 


Mg XI(i) 


9.23 


0.2±0.1 


Mg XI(f) 


9.31 


0.3 ±0.1 


Ne IX(r) 


13.45 


2.7 ±0.3 


Ne IX(i) 


13.55 


0.9 ±0.2 


Ne IX(f) 


13.70 


1.7 ±0.3 


Fe XIX 


13.52 


1.4 ±0.3 






0.68 confidence 


O VII(r) 


21.6 


< 1.4 


O VII(i) 
O VII(f) 


21.8 
22.1 


< 1.2 
<2.1 






0.90 confidence 


O VH(r) 


21.6 


<2.3 


O VII(i) 


21.8 


< 2.1 


O VII(f) 


22.1 


< 3.5 


Ne X Lya 


12.13 


1 3+" 4 

a '-0.6 
14 2+ 1 - 4 


FeXVn 
O VIII Lya 


15.01 
18.97 



NOTE. — Single line fluxes and He-like triplets obtained from the model 
of the MEG± 1 and HEG± 1 data using the combination of a bremsstrahlung 
model for the continuum and delta profiles for the lines. For the O VII triplet 
we include the flux estimate at 2 confidence levels; 68% and 90%. In the case 
of Ne IX we have also added a delta profile for Fe XIX which is blended with 
the Ne line. The fluxes are not corrected for the absorption column density. 



TABLE 7 

Line flux ratios and derived electron densities for the 
HDE 245059 binary 



0.68 confidence 0.90 confidence 

Triplet R n e (cm _J ) R n e (cnT 3 ) 

NelX 2.0±0.8 5±5 x 10 11 2.0±1.1 < 2 x 10 12 

MgXI 1.4 ±0.6 l +3 > 5 x 10 13 1.4± 1.0 l+* 8 xl0 13 

Si XIII 2.3 ±0.9 <5xl0 13 2.3 ±1.6 < 3 x 10 14 



NOTE. — Line flux ratios R = f/i and electron densities derived from 
the He-like triplets. Calculations were made for the grating spectrum of the 
binary (HE G and MEG first order) at two confidence levels, 68% and %9(). 
See section l4.4. ll for discussion. 
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Line fluxes for HDE 245059 binary components 



Kraus, A., Ireland, M., Martinache, E, Lloyd, J., & Hillenbrand, L. 2008, 

5'' ! Spitzer Conference. New Light on Young Stars: Spitzer's view of 

Circumstellar Disks. Pasadena, CA, USA. 
Lammer, H.,Selsis, F., Ribas, L, Guinan, E. F.. Bauer, S. J. & Weiss, W. W. 

2003, ApJ, 598, 121 
Lammer, H., Kulikov, Y. N. & Lichtenegger, H. I. M. 2006, Space Sci. Rev., 

122,189 

Lemen, J., R., Mewe, R., Schrijver, C, J. & Fludra, A. 1989, ApJ, 341, 474 
Li, J., Kastner, J.H., Prigozhin, G.Y, & Schulz, N. S. 2003, ApJ, 590, 586 
Li, J., Kastner, J.H., Prigozhin, G.Y, Schulz, N. S., Feigelson, E. D. & 

Getman, K. V. 2004, ApJ, 610, 1204 
Maddalena, R. J., Morris, M. 1987, ApJ, 323, 179 

McLean, I. S., Graham, J. R., Becklin, E. E., and Figer, D. F. 2000, Proc. 

SPIE Vol. 4008 
Murdin, P. & Penston, M.V. 1977, MNRAS, 181, 657 

Ness, J. U., Brickhouse, N. S., Drake. J. J. & Huenemoerder. D. P., 2003, 
ApJ, 598, 1277 

Neuhauser, R., Sterzik, M. E, Schmitt, J. H. M. M., Wichmann, R. & Krautter, 

J. 1995, A&A, 297, 391 
Oliveira, J. M., Jeffries, R. D., van Loon, J. T. & Rushton, M. T. 2006, 

MNRAS, 369, 272 
Padgett, D. L. 1996, ApJ, 471, 847 

Porquet, D.. Mewe, R., Dubau, J, Raassen, A. J. J., & Kaastra, J. S. 2001, 

yCat, 337,61113 
Rebull, L. M., Wolff, S. C, & Strom, S. E. 2004, AJ, 127, 1029 
Sacco, G. G., Franciosini, E., Randich, S. & Pallavicini, R. 2008, A&A, 488, 

167 

Schmitt, J. H. M. M., Robrade, J., Ness, J.U., Favata, F. & Stelzer, B. 2005, 

A&A, 432L, 35 
Siess, L., Dufour, E., and Forestini, M. 2000, A&A, 358, 593S 
Skinner, S. L., Brown, A. & Walter, F. M. 1991, AJ, 102, 1742 
Skinner, S. L. 1993, ApJ, 408, 660 

Skinner, S. L., Gudel, Koyama, K. & Yamaguchi, S. 1997, ApJ, 486, 886 
Smith, K., Audard, M., Gudel, M., Skinner, S. & Pallavicini, R. 2005, 

Proceedings 13th Cambridge Workshop on Cool Stars, Stellar Systems and 

the Sun 

Smith, D. S. & Scalo, J. 2007, Planet. Space Sci., 55, 517 
Stauffer, J. R. & Hartmann, L.W. 1986, PASP, 98, 1233 
Stelzer, B., Neuhauser, R. & Hambaryan, V. 2000, A&A, 356, 949 
Stelzer, B. & Schmitt, J. H. M. M. 2004, A&A, 418, 687 
Stone, R. C. & Taam, R. E. 1985, ApJ, 291, 183 

Telleschi, A., Gudel, M., Brigs, K. R., Audard, M., Ness, J. U. & Skinner, S. 

L. 2005, ApJ, 622, 653 
Telleschi, A., Gudel, M., Brigs, K. R., Audard, & M., Palla, F. 2007a, A&A, 

468, 425 

Telleschi, A., Gudel, M., Brigs, K. R., Audard, & M., Scelsi, L. 2007b, A&A, 
468, 443 

Testa, P., Drake, J. J. & Peres, G. 2004, ApJ, 617, 508 

Tsuboi, Y, Koyama, K, Murakami, H., Hayashi, M., Skinner, S. & Ueno, S. 

1998, ApJ, 503, 894 
Tsunemi, H., Mori, K, Miyata, E., Baluta, C, Burrows, D. N., Garmire, G. P., 

Chartas, G. 2001, ApJ, 554, 496 
Wizinowich, P. L., Acton, D. S., Lai, O., Gathright, J., et al. 2000, Proc. SPIE 

Vol. 4007 







Flux (10" 


? T 

photons cm s ) 


Ion 


A (A) 


north 


south 


MgXII 


8.42 


0.5 ±0.1 


0.1 ±0.1 


Ne X Lya 


12.13 


3.7 ±0.4 


2.8 ±0.4 


Ne X Ly/3 


10.24 


0.5±0.1 


0.2±0.1 


O VIII Lya 


18.97 


7.4+f 


5.7 + f 



NOTE. — Single line flux obtained for each star of the binary from the 
model of the MEG± 1 data using the combination of a bremsstrahlung model 
for the continuum and a delta function for the line profiles. The fluxes are not 
corrected for the absorption column density. 



